Amendments of the Claims: 



A detailed listing of all claims in the application is presented below. This listing of claims 
will replace all prior versions, and listings, of claims in the application. All claims being 
currently amended are submitted with markings to indicate the changes that have been made 
relative to immediate prior version of the claims. The changes in any amended claim are being 
shown by strikethrough (for deleted matter) or underlined (for added matter). 

13. (Cancelled) 

14. (Cancelled) 

15. (Currently Amended) A practical numerical computer-implemented method for reliably 

computing a dynamical decomposition point of a stable equilibrium point for large-scale 
nonlinear systems, comprising the steps of: 

a) given a stable equilibrium point x s . 



b) moving along a search path cp t (x s ) = {jc, + t x i, t e 5R + } starting from x s and 

detecting an exit point, x a , at which said search path <p t (x s ) ex * ts a 
stability boundary of a stable equilibrium point x s ; 

c) using said exit point x^ as an initial condition and integrating a nonlinear 

system (4.2) to an equilibrium point x d ; and 

d) computing said dynamical decomposition point with respect to the stable 

equilibrium point x St wherein said search path dir e ction S ~ is e x d .; and 

e) displaying the dynamical decomposition point. 
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16. (Previously Presented) The method of claim 15, wherein a method for computing said exit 

point of the nonlinear system (4.3) comprises the step of moving along said search path 

9t ( x s) = { x s + t x s> r e 9? + } starting from x s and detecting said exit point x^ , which is 
a first local maximum of an objective function C(x) along said search path <p t (x s ) . 

17. (Previously Presented) The method of claim 15, wherein a method for computing a 

dynamical decomposition point comprises the steps of: 

a) using said exit point x a as an initial condition and integrating a nonlinear 

system (4.2) to a first, local minimum of a norm ||F(x)|| along the 
corresponding trajectory, where F(x) is a vector field of (4.2), and letting 
the point at which the first local minimum of ||F(x)|| occurs be denoted 
jc° , and is called the minimum distance point (MDP); and 

b) using said MDP x] as an initial guess and solving a set of nonlinear algebraic 

equations of said vector field (4.2) F(x) = 0 , wherein a solution is x d , and 
a dynamical decomposition point with respect to the local optimal solution 
x s and said search path <p t ( x s ) * s x d • 

18. (Previously Presented) The method of claim 15, wherein a method for computing said exit 

point with respect to a stable equilibrium point of the nonlinear system (4.2) and a search 
vector comprises the step of computing an inner-product of said search vector and the 
vector field of system (4.3) at each time step, by moving along said search path 

<p t (x s )= {x s + tx s, t € W + } starting from x s and at each time-step, computing an inner- 
product of said search vector s and vector field F(x) , such that when a sign of said 
inner-product changes from positive to negative, said exit point is detected. 



19. (Previously Presented) The method of claim 15, wherein a method for computing said exit 

point of nonlinear system (4.3) with respect to a stable equilibrium point and a search 
vector comprises the steps of: 

a) moving from said stable equilibrium point along said search vector until an 

inner-product of said search vector and the vector field of system (4.3) 
changes sign between an interval [t\, t{\\ 

b) applying a linear interpolation to an interval [t\ 9 feL which produces an 

intermediate time to where an interpolated inner-product is expected to be 
zero; 

c) computing an exact inner-product at /o> such that if said value is smaller than a 

threshold value, said exit point is obtained; and 

d) if said inner-product is positive, then replacing t\ with to, and otherwise 

replacing h with to and going to step b). 

20. (Previously Presented) The method of claim 15, wherein a method for computing a minimum 

distance point (MDP) of the nonlinear system (4.2) satisfying conditions (CI) and (C2) 
comprises the steps of: 

a) using said exit point as an initial condition and integrating the nonlinear system 

(4.2) for a few time-steps, and letting the end point be denoted as the 
current exit point; 

b) checking convergence criterion, and, if a norm of said current exit point 

obtained in step a) is smaller than a threshold value, then declaring said 
point as said MDP and stopping the process, otherwise, going to step c); 
and 

c) drawing a ray connecting a current exit point on a trajectory and a local optimal 

solution (equivalently, a stable equilibrium point), replacing said current 
exit point with a corrected exit point, which is a first local maximal point 
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of objective function along said ray, starting from the stable equilibrium 
point, and assigning this point to said exit point and going to step a). 



21. (Previously Presented) The method of claim 15, wherein a method for computing said 

dynamical decomposition point of the nonlinear system (4.2) satisfying conditions (CI) 
and (C2) with respect to a stable equilibrium point x s and a search vector $ , comprises 
the steps of: 



a) moving along said search path <p , (x s ) = \x s + t x s ,t e 9* + } starting from x s and 

detecting a moment that an inner-product of said search vector s and the 
vector field F(x) of nonlinear system (4.2) changes sign, between an 

interval j/ 1? f 2 ] , stopping this step if t x is greater than a threshold value and 
reporting that there is no adjacent local optimal solution along this search 
path, otherwise, going to step b); 

b) applying linear interpolation to said interval [fj ,f 2 ] , which produces an 

intermediate time t 0 where said interpolated inner-product is expected to 
be zero, computing an exact inner-product at t 0 , and if said value is 
smaller than a threshold value, said exit point is obtained, and going to 
step d); 

c) if said inner-product is positive, then replacing t x with t 0 , and otherwise 

replacing t 2 with t 0 and going to step b); 

d) using said exit point as an initial condition and integrating a nonlinear system 

(4.2) for a few time-steps, and letting the end point be denoted as the 
current exit point; 

e) checking convergence criterion, and if a norm of said point obtained in step d) 

is smaller than a threshold value, then declaring said point as the MDP and 
going to step g), otherwise going to step e); 



7 



f) drawing a ray connecting a current exit point on said trajectory and a stable 

equilibrium point, replacing said current exit point with a corrected exit 
point which is a first local maximal point of objective function along said 
ray starting from the stable equilibrium point, and assigning this point to 
said exit point and going to step d); and 

g) using said MDP as an initial guess and solving a set of nonlinear algebraic 

equations of the vector field (4.2) F(x) = 0 , wherein a solution is t d , such 
that said DDP with respect to a stable equilibrium point x s and search 
vector s is x d . 

22. (Cancelled) 

23. (Previously Presented) A hybrid search computer-implemented method for obtaining a local 

optimal solution of general unconstrained nonlinear programming problem (4.1) starting 
from any initial point, comprising the steps of: 

a) given an initial point x 0 ; 

b) integrating a nonlinear dynamical system described by (4.2) that satisfies 

conditions (CI) and (C2) from said initial point x Q to obtain a trajectory 

0,(x o ) for n time-steps, n being an integer, and recording the last time-step 

point as the end point, and if it converges to a local optimal solution, then 
stopping, otherwise, going to step (c); 

c) monitoring a desired convergence performance criterion along the trajectory 

<f> t {x 0 ) in terms of the rate of decreasing values in the objective function 
under study, and if the desired criterion is satisfied, then using the end 
point of trajectory ^(x 0 ) as the initial point and going to step (b), 
otherwise, going to step (d); and 
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d) applying an effective local optimizer (i.e. 9 a method to find a local optimal 
solution) from said end point in step (b) to continue the search process, 
and if it finds a local optimal solution, then displaying the local optimal 

solution and stopping, otherwise, setting the end point of trajectory 0,(x o ) 

as the initial point, namely x 0 , and going to step (b). 

24. (Cancelled) 

25. (Previously Presented) The method of claim 44-23, wherein in step (a) the initial point x 0 is 

given by a method comprising the steps of: 

a) moving along a search path starting from a local optimal solution x opt and 

applying a DDP search method to compute a corresponding DDP, and if a 
DDP can be found, then going to step (b), otherwise, trying another search 
path and repeating this step; 

b) letting said DDP be denoted as x d , and if x d has previously been found, then 

going to step a), otherwise going to step c); and 

c) setting x 0 = x opt + (1 + € ){x d - x opt ) where e is a small number, and applying a 

hybrid search method starting from x 0 to find a corresponding adjacent 
local optimal solution. 

26. (Previously Presented) A computer-implemented method for obtaining a global optimal 

solution of unconstrained nonlinear optimization problems, comprising the steps of: 

a) choosing a starting point; 

b) applying the hybrid search method of claim 23 using said starting point to find 

a local optimal solution x® ; 

c) setting; = 0, V s = {x° s } , Vj ew = {x° s }. and V d = {#} ; 
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d) wherein setVjJ^ = {<j>) and for each local optimal solution in V n J ew (i.e., x J s ), 

performing steps (e) through (k); 

e) defining a set of search vectors Sf , / = 1,2, . . my, and setting i = 1; 

f) if / > ntj , then going to step (1); otherwise, applying a DDP search method 

along the search vector Sf to find a corresponding dynamical 
decomposition point (DDP), and if a DDP is found, then going to step (g), 
otherwise, setting i = i + 1 and going to step (f); 

g) letting the found DDP be denoted as , , checking whether it belongs to the 

set V d , i.e., x J d , e V d ?, and if it does, then setting i = i + 1 and going to step 
(f), otherwise setting V d = V d u and going to step (h); 

h) for the DDP x J dJ , performing steps (i) through (j) to find a corresponding 

adjacent local optimal solution; 

i) letting x J Ql = x J s + (l + e)(x J di - x J s ) , where s is a small number; 

j) applying said hybrid search method using x J 0 i as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as x' sj ; 

k) checking whether x' sj has been found before, i.e., x' sj e V s ?, and if it has 

already been found, then setting / = / + 1 and going to step (f), otherwise, 
setting V s = V s u } and = u {<,.} and setting / = / + 1 and 
going to step (f); 

1) examining the set of all newly computed local optimal solutions, V n J e + J , and if 
V+£ is non-empty, then setting j= j + 1 and proceeding to step (d), 
otherwise proceeding to the next step; and 
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m) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their corresponding objective function values in set V s 
and; 

n) displaying the global optimal solution. 

27. (Cancelled) 

28. (Currently Amended) A computer-implemented method for obtaining a global optimal 

solution of a_constrained nonlinear programming problem (4.5), comprising two phas e s 
the steps of : 

a) Phase I: finding all feasible components of the constrained nonlinear 

programming p roblem; and 

b) Phase II: finding all local optimal solutions for the constrained nonlinear 

programming problem in each feasible component; 

c) chosing choosing a global optimal solution for the constrained nonlinear 

programming problem from the local optimal solutions found in step (b); 
and 

d) displaying the global optimal solution for the constrained nonlinear 

programming problem . 

29. (Previously Presented) The method of claim 28, comprising the steps of: 

a) approaching a path-connected feasible component of a constrained 

optimization problem (4.5); and 

b) escaping from said path-connected feasible component and approaching 

another path-connected feasible component of said constrained 
optimization problem (4.5). 

30. (Previously Presented) The method of claim 28, comprising the steps of: 
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a) approaching a stable equilibrium manifold of a nonlinear dynamical system 

(4.9) satisfying conditions (CI -1) and (CI -2); and 

b) escaping from said stable equilibrium manifold and approaching another stable 

equilibrium manifold of said nonlinear dynamical system (4.9) satisfying 
conditions (C 1-1) and (CI -2). 

31. (Previously Presented) The method of claim 28, comprising the steps of: 

a) starting from a point in a feasible component and approaching a local optimal 

solution located in said feasible component of an optimization problem 
(4.7); and 

b) escaping from said local optimal solution and approaching another local 

optimal solution of said feasible component of said optimization problem 
(4.7). 

32. (Previously Presented) The method of claim 28, comprising the steps of: 

a) in a deterministic manner, first finding all stable equilibrium manifolds of a 

nonlinear dynamical system that satisfies conditions (Cl-1) and (CI -2); 

b) in a deterministic manner, finding all stable equilibrium points of a nonlinear 

dynamical system that satisfies conditions (C2-1) and (C2-2); and 

c) then from said stable equilibrium point, finding a global optimal solution. 

33. (Previously Presented) The method of claim 32, comprising the steps of: 

a) given a feasible solution of constrained nonlinear programming problem (4.5); 

b) finding a stable equilibrium point of a nonlinear dynamical system (4.10) 

satisfying conditions (C2-1) and (C2-2); 

c) moving from said stable equilibrium point to a dynamical decomposition point, 

in order to escape from a local optimal solution; and 
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d) approaching another stable equilibrium point of said nonlinear dynamical 

system (4.10) satisfying conditions (C2-1) and (C2-2) in the same path- 
connected feasible component, via said dynamical decomposition point. 



34. (Cancelled) 

35. (Cancelled) 

36. (Previously Presented) A hybrid search method for Phase I of claim 28, wherein one 

effective local method is combined with said dynamical trajectory method, comprising 
the following steps to find a feasible solution of the constrained optimization problem 



a) given an initial point x 0 ; 

b) integrating a nonlinear dynamical system described by (4.9) that satisfies 

conditions (Cl-1) and (CI -2) from said initial point to obtain a trajectory 

$,{x 0 ) for n time-steps, n is an integer, and recording the last time-step 
point as the end point, and if it converges to a feasible solution, then 
stopping, otherwise, going to step (c); 

c) monitoring a desired convergence performance criterion in the objective 

function ||#(*)||, a vector norm of M*) in (4.7), and if the desired 



criterion is satisfied, then using the end point of trajectory $ t [x 0 ) as the 
initial point and going to step (b), otherwise, going to step (d); and 

d) applying an effective method from the said end point in step (b) to continue the 
search process, and if it finds a feasible solution of (4.7) then stopping, 

otherwise, setting the end point of trajectory 0, (jc 0 ) in step (b) as the initial 

point and going to step (b). 



(4.7); 




13 



37. (Previously Presented) A hybrid search method for Phase II of claim 28, comprising the 

following steps to find a local optimal solution of the constrained optimization problem 
(4.7). 

a) given a feasible point x 0 ; 

b) integrating a nonlinear dynamical system described by (4.10) that satisfies 

conditions (C2-1) and (C2-2) from said initial point to obtain a trajectory 

<j> t (x Q ) for n time-steps, n is an integer, and recording the last time-step 
point as the end point, and if it converges to a local optimal solution, then 
stopping, otherwise, going to step (c); 

c) monitoring a desired convergence performance criterion in an objective 

function C(x) (©) of (4.5)., and if the desired criterion is satisfied, then 

using the end point of trajectory ^(a; 0 ) as the initial point and going to 
step (b), otherwise, going to step (d); and 

d) applying an effective method from the said end point in step (b) to continue the 

search process, and if it finds a feasible solution of (4.7) then stopping, 

otherwise, setting the end point of trajectory ^(x 0 ) in step (b) as the initial 
point and going to step (b). 

38. (Previously Presented) The method of claim 30, comprising the steps of: 

a) choose a starting point; 

b) initialization j = 0; 

c) applying a hybrid search method for Phase I using said starting point to find a 

feasible point in a (path-connected) feasible component, setting a point so 
found as an initial point, and applying the hybrid search method for Phase 
II to find a local optimal solution x J sfi ; 
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d) starting from said initial point x J s Q , applying a numerical method for Phase II to 

find all local optimal solutions in said feasible component and recording 
them as the set V s j and set V s = V s u Vj ; 

e) for the local optimal solution x/ 0 , defining a set of search vectors 

8/ J - l>2,...,£y ; 

f) if i > Aj, then going to step (1), otherwise, going to step (g); 

g) for each search vector S{ , applying a reverse-time trajectory method to find a 

point lying in an unstable equilibrium manifold of a system (4.9), and if a 
point is found, letting it be denoted as x l dJ and going to step (h), otherwise 
setting i' = i + 1 and going to step (f); 

h) setting 4 y = x j s + (/ + s)(x dj - x J s ), where e is a small number and x J s is a local 

optimal solution selected in step (e), and applying the hybrid search 
method for Phase II using x* 0 d as the initial point to find a point lying in a 
stable equilibrium manifold of system (4.9) satisfying conditions (Cl-1) 
and (CI -2), and letting the solution be denoted as x° sj ; 

i) starting from said initial point x° sj and applying the hybrid search method of 

step (b) of claim 28 to find a local optimal solution in said (path- 
connected) feasible component x J s *J ; 

j) checking whether x J s ^ has been found before, i.e. x J s * 0 x € V s ?, and if it has been 

bound before (Le. 9 the said feasible component has been visited before), 
then setting i = / + 1 and going to step (f), otherwise setting 

K = V s u {x J St0 } and V£ = V£ u [x$ ) and going to step (k); 
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k) starting from said initial point x J s * 0 x , applying a numerical method for Phase II 
to find all local optimal solutions in said feasible component and recording 
them in the set F/ +1 , set V s = V s u V s J+] and setting i = i + 1 and going to 
step (f); 

1) examining the set of all newly computed local optimal solutions , and if 
Vjjj is empty, then going to the next step, otherwise, setting^ =y + 1 and 
going to step (e); and 

m) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their objective function values. 

39. (Previously Presented) A numerical method for Phase II of claim 28, comprising the 

following steps to find all the local optimal solutions in said path-connected feasible 
component: 

a) given an initial point lying in said path-connected feasible component; 

b) initialization; 

c) applying a hybrid search method for Phase II using the initial point to find a 

local optimal solution x° s , and setting j = 0, V s = n = [x° s ] and 

d) for each local optimal solution in Vj^ (i.e., x j s \ performing steps (e) through 

(k); 

e) defining a set of search vectors Sf 9 i = 1,2,..., m j , and setting i = 1; 

f) if i > m h then going to step (1), otherwise, applying a DDP search method, 

wherein a nonlinear dynamical system (4.10) satisfies conditions (C2-1) 
and (C2-2) to find a corresponding dynamical decomposition point (DDP), 
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and if a DDP is found, then going to step (g), otherwise, setting i = i + 1 
and going to step (f); 

g) letting the found DDP be denoted as x J dJ , checking whether it belongs to the 

set V d , i.e., , e V d ?, and if it does, then setting / = / + 1 and going to step 

(f), otherwise setting V d = V d u jx^. J and going to step (h); 

h) for the DDP x J di , performing steps (i) through (j) to find a corresponding 

adjacent local optimal solution; 

i) setting x J 0i = x J s + (l + s){x J di - x J s ^ , where s is a small number; 

j) applying a hybrid search method using x J 0 i as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as x l s j ; 

k) checking whether x' sj has been found before, i.e. x l s j € V s ?, and if it has 

already been found, then setting / = i + 1 and going to step (f), otherwise, 
setting V s ~V s u {x' tJ } and = u {x l sj ) and setting / = i + 1 and 
going to step (f); and 

1) examining the set of all newly computed local optimal solutions Vj^ u , and if 
Kw,u * s non-empty, then settings =7 +1 and proceeding to step (d), 
otherwise outputting all the local optimal solutions in said path-connected 
feasible component contained in set V s and stopping the process. 

40. (Previously Presented) A computer-implemented method for obtaining the global optimal 
solution of constrained nonlinear programming problems, comprising the steps of: 

a) using a transformation technique to transform a constrained optimization 

problem into an unconstrained optimization problem, then applying the 
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following steps to find the global optimal solution of the unconstrained 
optimization problem; 

b) choosing a starting point; 

c) apply a hybrid search method using said starting point to find a local optimal 

solution jc° ; 

d) setting^ = 0, V s = {x° s } , Vj ew = {x° s } and V d = {«>} ; 

e) setting = {</>} and for each local optimal solution in (i.e., x J s ), 

performing steps (e) through (k); 

f) defining a set of search vectors Sf , i = 1 ,2, . . ny, and setting i = 1 ; 

g) if / > rrij , then going to step (1), otherwise, applying the DDP search method of 

claim 15 along the search vector Sf to find a corresponding dynamical 
decomposition point (DDP), and if a DDP is found, then going to step (g), 
otherwise, setting / = / + 1 and going to step (f); 

h) letting the found DDP be denoted as x J dJ , checking whether it belongs to the 

set V d , i.e., x J dJ e V d ?, and if it does, then setting / = / + 1 and going to step 
(f), otherwise setting V d = V d u jxj , j and going to step (h); 

i) for the DDP x J d i , performing steps (i) through (j) to find a corresponding 

adjacent local optimal solution; 

j) setting x J Qi = x j s + (l + s){x J di - x J s } , where s is a small number; 

k) applying a hybrid search method using jc^. as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as x l s j ; 
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1) checking whether x' sj has been found before, i.e. x' sj e V s ?, and if it has 

already been found, then setting i = / + 1 and going to step (f), otherwise, 
setting V s = V s u {<,} and V£ l = Vj e + J u {<,} and setting i = / + 1 and 
going to step (f); 

m) examining the set of all newly computed local optimal solutions, , and if 
Vj e l l is non-empty, then setting j= j+l and proceeding to step (d), 
otherwise proceeding to the next step; aad 

n) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their corresponding objective function values in set V s ; 
and 

o) displaying the global optimal solution. 
Please cancel claim 41 in its entirety. 
41. (Cancelled) 
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